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Abstract. Cooling of neutron stars (NSs) with the cores 
composed of neutrons, protons, and electrons is analyzed. 
The main cooling regulators are discussed: opening of di- 
rect Urea process in a NS central kernel; superfluidity of 
nucleons in NS interiors; surface layers of light (accreted) 
elements; strong surface magnetic fields. An emphasis is 
paid on the cooling scenario with strong 1 Sq pairing of pro- 
tons and weak 3 P2 pairing of neutrons in the NS core, as 
well as strong x So pairing of neutrons in the NS crust. The 
theory predicts three types of isolated cooling middle-aged 
NSs with distinctly different properties: low-mass, slowly 
cooling NSs; medium-mass, moderately cooling NSs; mas- 
sive, rapidly cooling NSs. The theory is compared with 
observations of isolated NSs — pulsars and radio quiet 
NSs in supernova remnants. The constraints on physical 
properties of NSs which can be inferred from such a com- 
parison are outlined. 



1. Introduction 

Cooling of neutron stars (NSs) depends on the properties 
of dense matter in their crusts and cores. These properties 
are still poorly known and cannot be described unambigu- 
ously by contemporary theories. For instance, calculations 
of the equation of state (EOS) of the NS cores (e.g., Lat- 
timer & Prakash |2001 ) or the superfluid pro pertie s of NS 
cores and crusts (e.g., Lombardo & Schulze 2001 ) give a 
large scatter of results which depend on a model of strong 
interaction and a many-body theory employed. However, 
these properties can be studied by comparing the results of 
simulations of NS cooling with the observations of thermal 
emission from isolated NSs. We will describe some recent 
results of such studies. The history of NS c ooling theory 
is reviewed, for instance, by Yakovlev et al. ( 1999| ) 

For simplicity, we consider the NS models with the 
cores composed mainly of neutrons with an admixture of 
protons and electrons. We discuss the basic ideas of the 
cooling theory, and the main cooling regulators, paying 
special attention on the effects of superfluidity of nucleons 
in NS interiors. Our consideration will be based on recent 
simulations of the cooling of superfluid NSs by Kaminker 



et al. fl200l| ), Potekhin & Yako vlev (p00l| ), Yakovlev et al. 
( |2001b| ), and Kaminker et al. ( |2002| , hereafter KYG). 



2. Equations of thermal evolution 

Neutron stars are born very hot in supernova explosions, 
with internal temperature T ~ 10 11 K, but gradually cool 
down. In about 30 s after the birth a star becomes fully 
transparent for neutrinos generated in its interiors. We 
consider the cooling in the following neutrino-transparent 
stage. The cooling is realized via two channels, by neu- 
trino emission from the entire stellar body and by heat 
transport from the internal layers to the surface result- 
ing in the thermal emission of photons. For simplicity, we 
neglect the possible reheating mechanisms (frictional dis- 
sipation of rotational energy, or Ohmic decay of internal 
magnetic field, or the dissipation associated with weak de- 
viations from the chemical equilibrium as reviewed, e.g., 
by Page |l998|) . 

The internal structure of NSs can be regarded 
as te mperature-independent (e.g., Shapiro & Teukolsky 
1983). The general relativistic equations of thermal evo- 



lution include the energy and flux equations obtained by 
Thorne (1977). For a spherically symmetric NS, 
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where Q is the neutrino emissivity [erg cm~ 3 s _1 ], c v is the 
heat capacity per cm 3 [erg cm -3 K _1 ], k is the thermal 
conductivity [erg cm -1 s _1 K _1 ], and L r is the "local lu- 
minosity" [erg s _1 ] defined as the non- neutrino heat flux 
transported through a sphere of radius r. The gravita- 
tional mass m(r) and the metric function $(r) are deter- 
mined by the stellar model. 



It is conventional (e.g., Gudmundsson et al. 1983) to 



subdivide the calculation of heat transport in the neutron- 
star interior (r < i?b) and in the outer heat-blanketing 
envelope (i?b < t < R), where R is the stellar radius, 
and the boundary radius Ri, corresponds to the density 
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Pb ~ 10 10 g cm -3 (~ 100 meters under the surface). In 
the presence of superstrong surface magnetic field, it is 
reasonable to shift ph to the neutron drip density, 4 x 



Table 1. 

KYG) 



NS models employing EOSs A and B (from 



10 11 g cm" 3 (Potekhin & Yakovlev 2001). The thermal 



structure of the blanketing envelope is studied separately 
in the stationary, plane-parallel approximation to relate 
the effective surface temperature T s to the temperature Xb 
at the inner boundary of the envelope. We discuss the T s - 
Tb relation in Sect. 3. It is used as the boundary condition 
for solving Eqs. (0) and (||) at r < i?b- 

The effective temperature determines the photon lu- 
minosity L 1 — 47rcri? 2 T s 4 (t). Both, L 1 and T s refer to 
the locally-flat reference frame on the surface. A dis- 
tant observer would register the "apparent" luminosity 
L^ 3 = L 7 (l — r g /R) and the "apparent" effective temper- 
ature T s °° = T s y/1 - r s /R, where r g = 2GM/c 2 is the 
Schwarzschild radius and M = m{R) is the total gravita- 
tional mass. 

The main goal of the cooling theory is to calculate 
cooling curves, T£°(t), to be compared with observations. 
One can distinguish three main cooling stages: (i) the in- 
ternal relaxation stage (t <, 10-100 yr; Lattimer et al. 



1994; Gnedin et al. 2001), (ii) the neutrino cooling stage 
(the neutrino luminosity L v 3> L<y, t ^ 10 5 yr), and (iii) 
the photon cooling stage (L v <C L 7 , 10 5 yr). 

After the thermal relaxation, the redshifted temper- 
ature T(t) = T(r,t) e*( r ) becomes constant throughout 
the stellar interior. Then Eqs. ([!]) and (||) reduce to the 
equation of global thermal balance, 



C(T) 



dT 



dt 

L~(f) = 47T 

C(f) = 4tt 



L~(T)-i~(T a ), 

Rb drr 2 Q(f) e 2 " 6 
o y/l~2Gm/(W^)' 
R " drr 2 c v (T) 
o y/1 - 2Gm/(cV]' 



(3) 
(4) 

(5) 



where C is the total NS heat capacity, and is the total 
neutrino luminosity. 

3. Physics input 

We describe the results obtained using a full y rela tivis- 
tic non- isothermal cooling code (Gnedin et al. [2001 ). The 
code solves the heat diffusion equations (Sect. 2) in the 
NS interior. 

The physics unput is as follows. The EOS in the NS 
crust is taken from Negele & Vautherin ( 1973 ) (assuming 
spherical atomic nuclei everywhere in the crust). The core- 
crust boundary is placed at the density 1.5 x 10 14 g cm~ 3 . 
In the core, we use two phenomenologica l EO Ss (EOS A 
and EOS B) proposed by Prakash et al. ( |l988|) . 

EOS A is model I of Prakash et al. (|1988|) with the 



Model 


Main parameters 


EOS A 


EOS B 


Maximum 


M max /M 


1.977 


1.73 


mass 


Pcmax/10 14 g Cm" 3 


25.75 


32.5 


model 


R km 


10.754 


9.71 


Direct Urea 


A/d/Mq 


1.358 


1.44 


threshold 


po/10 14 g cm" 3 


7.851 


12.98 


model 


7? km 


12.98 


11.54 



compression modulus of saturated nuclear matter K = 
240 MeV. EOS B corresponds to K = 180 MeV and to 



the simplified form of the symmetry energy suggested by 
Page & Applegate ( |1992 ). It has been used in many papers 
(e.g., Page & Applegate 1992| , Yakovlev et al. |199S| , pOOla , 
and references therein). 

The masses, central densities, and radii of two stellar 
configurations for EOSs A and B are given in Table 1. 
The first configuration corresponds to the maximum-mass 
NS. The values of M max reveal that EOS A is stiff, and 
EOS B is moderate. The second configuration has a central 
density p c at which the direct Urea process switches on 
(the process is allowed at p c > pd, which corresponds to 
M > Md). EOS B is based on a smaller symmetry energy 
and opens the direct Urea process at a higher p. 

The cooling code includes all important neutrino emis- 
sion processes in the NS core (direct and modified Urea 
processes, neutrino bremsstrahlung in nucleon-nucleon 
scattering, neutrino emission due to Cooper pairing of 
nucleons) and in the crust (plasmon decay, neutrino 
bremsstrahlung due to scattering of electrons off atomic 
nuclei, electron-positron annihilation into neutrino pairs, 
neutrino emission due to Cooper pairing of neutrons in the 
inner crust). The se proc esses are reviewed, for instance, by 
Yakovlev et al. ( 2001a ). The effective masses of protons 
and neutrons in the core and free neutrons in the crust 
are taken to be 0.7 of the bare nucleon masses. The values 
of the thermal conductivity in the NS crust and core are 
the same as used by Gnedin et al. (|200l| ). 

The T s -Tb relationship is taken from Potekhin et al. 
( |1997| ) and Potekhin & Yakovlev ( [200 1|) . It is relevant 
either for the surface layers made of iron (without any 
magnetic field and with the dipole surface magnetic fields 
B ^ 10 15 G), or for the non-magnetic surface layers con- 
taining light elements. In addition, we have modified the 
relationship by incorporating our new results which take 
into account the combined effects of accreting envelopes 
and magnetic fields. It is assumed that the surface mag- 
netic field induces an anisotropic heat transport in the 
heat-blanketing envelope but does not violate the isotropic 
(radial) heat diffusion in the deeper NS layers. This pre- 
serves the scheme of solving the radial heat diffusion equa- 
tions (Sect. 2). In this case, T s varies over the NS surface, 
and the cooling theory deals with the mean surface tem- 
perature T s that determines the NS photon luminosity. 
For brevity, we denote T s as T s . It is also assumed that 
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Table 2. Parameters of superfluid models in Eq. 
KYG 



from 



Pair- 


Mo- 


To/10 9 K 


kg 


ing 


del 




fm 


'S 


lp 


20.29 





'So 


2p 


17 





% 


3p 


14.5 





'So 


Ins 


10.2 





'So 


2ns 


7.9 





'So 


3ns 


1800 





3 P 2 


lnt 


6.461 


1 



fm" 1 



k 2 

fm" 1 



fc 3 

fnT 



1.117 

1.117 

1.117 

0.6 

0.3 

21 

1.961 



1.241 

1.329 

1.518 

1.45 

1.45 

1.45 

2.755 



0.1473 

0.1179 

0.1179 

0.1 

0.01 

0.4125 

1.3 



a NS may have a hydrogen atmosphere even if the heat- 
blanketing envelope is mostly made of iron. The majority 
of cooling curves are calculated for the non-magnetized 
heat-blanketing envelopes made of iron. The exceptions 
are discussed in Sect. 6.3. 

We will focus on the effects of nucleon superfluidity on 
NS cooling. The superfluid properties of NS matter are 
characterized by the density-dependent critical tempera- 
tures T c (p) of nucleons. Microscopic theories predict (e.g., 
Lombardo & Schulze 2001, Yakovlev et al. [1999 ) super- 
fluidities of three types: singlet-state ( 1 So) pairing of neu- 
trons (T c — T cns ) in the inner crust and outermost core; 
1 So proton pairing (T c — T cp ) in the core; and triplet-state 
( 3 P2) neutron pairing (T c — T cnt ) in the core. Superflu- 
idity of nucleons suppresses neutrino processes involving 
nucleons, initiates a specific mechanism of neutrino emis- 
sion a ssocia ted with Cooper pairing of nucleons (Flowers 
et al. 1976 ), and affects the nucleon heat capacity. All 
these effects are incorporated into the cooling code as de- 
scribed by Yakovlev et al. (1999, 2001a) and Gnedin et al. 
( |200lD . 

The microscopic calculations of superfluid critical tem- 
peratures T c (p) give a large scatter of results (e.g., Lom- 
bardo & Schulze 200 1[ ), but some common features are 
clear. For instance, T c (p) increases with p at sufficiently 
low densities (due to an increasing strength of the attrac- 
tive part of nucleon- nucleon interaction), reaches maxi- 
mum and then decreases (due to a short-range nucleon- 
nucleon repulsion), vanishing at a rather high density. 
For T cns (p), the maximum occurs at sub- nuclear densi- 
ties, while the switch off takes place at p ~ po, where 
po « 2.8 x 10 14 g cm -3 is the saturated nuclear matter 
density. For T cp (p) and T cnt (p), the maxima take place at 
a few po and the fall occurs at the densities several times 
higher. The maximum values of T c range from about 10 8 
K (or even lower) to (2 — 3) x 10 10 K, depending on the mi- 
croscopic theoretical model. The maximum values of T cnt 
are typically lower than those of T cp and T cns , due to the 
weaker nucleon-nucleon attraction in the 3 P2 state. 

Instead of studying T c as a function of p, it is often con- 
venient to consider T c as a function of the nucleon Fermi 
wavenumber k = k-p^ = (37t 2 7t.n) 1 ^ 3 , where tin is the num- 
ber density of nucleon species N=n or p (neutrons or pro- 
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lg p (g cm" 3 ) 

Fig. 1. Density dependence (from KYG) of the critical 
temperatures for three models lp, 2p, and 3p of the pro- 
ton superfluidity (dots-and-dashes) in the core (with EOS 
A); three models Ins, 2ns, and 3ns of neutron super- 
fluidity (solid, short-dashed, and long-dashed lines); and 
one model lnt of P2 neutron superfluidity (dots) used 
in cooling simulations. Vertical dotted lines indicate the 
neutron drip point, core-crust interface, and direct Urea 
threshold. 

tons). Moreover, instead of T c (k) one often considers A(fc), 
the zero-temperature superfluid gap. For the 1 So pair- 
ing, assuming BCS theory, one has A(fc) = T c (fc)/0.5669. 
In the case of 3 P2 neutron pairing, the gap depends on 
the orientation of nucleon momenta with respect to the 
quantization axis. Following the majority of papers, we 
adopt the 3 P2 pairing with zero projection of the to- 
tal angular momentum on the quantization axis. In that 
case A nt (fc) = T cnt (/c)/0.8416, where A nt (fc) is the mini- 
mum value of the gap on the neutron Fermi surface (e.g., 
Yakovlev et al. 1999| ). The dependence of the gaps on k 
for the selected superfluidity models is plotted in KYG. 

Taking into account large scatter of theoretical values 
of T c (k) it is instructive to consider T c (k) as unknown 
functions and analyze (Sect. 6) which of them are consis- 
tent with observations. In this paper, following KYG, we 
adopt the parameterization of T c of the form: 



T c = To 



(k - k f 



(k - k 2 f 



(6) 



(k-ko^ + kf (fc-fc 2 ) 2 + fc| ' 

for ko < k < k 2 ; and T c = 0, for k < kg or k > k 2 . The 
factor To regulates the amplitude of T c , fco and k 2 deter- 
mine positions of the low- and high-density cutoffs, while 
ki and k 3 specify the shape of T c (p). All wave numbers, 
k, fco, . . . &3 are expressed in fen . KYG verified that by 
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Fig. 2. Density dependence of the critical temperatures of 
superfluidity of protons (models lp, 2p, and 3p — solid, 
dashed, and dot-and-dashed lines) and neutrons (model 
lnt — dotted line) in the NS core for EOS A on a larger 
scale than in Fig. [l[ Vertical dotted line indicates direct 
Urea threshold. Horizontal line is the temperature of mat- 
ter, T = 3 x 10 8 K, adopted in Fig. |. 



tuning To, ko, ... A3, the parameterization accurately de- 
scribes numerous results of microscopic calculations. 

Following KYG we use three models of 1 So proton su- 
perfluidity, three models of neutron superfluidity, and 
one model of 3 P2 neutron superfluidity. The parameters 
of the models are given in Table 2, and the appropriate 
T c (p) are plotted in Figs. [I], ||, and [l| One has fc = for 
1 So pairing. At any given p it is reasonable to choose the 
neutron superfluidity ^So or 3 P2) with higher T c . 

Models Ins and 2ns of neutron pairing in the crust 
correspond to about the same, rather strong superfluidity 
(with maximum T cns w 7 x 10 9 K). Model 2ns has flat- 
ter maximum and sharper decreasing slopes in the wings 
(near the crust-core interface and the neutron drip point). 
Model 3ns represents a much weaker superfluidity, with 
maximum T cns « 2.4 x 10 9 K and a narrower density pro- 
file. 

The proton superfluidity curves lp, 2p, and 3p in Figs, 
[l] and are similar. The maximum values of T cp are about 
7 x 10 K for all three models. The models differ by the po- 
sitions of the maximum and decreasing slopes of T cp (p). 
The decreasing slope of model lp is slightly above the 
threshold density of the direct Urea process (for EOS A), 
while the slopes for models 2p and 3p are shifted to higher 
p. Our models lp, 2p, and 3p are typical for those micro- 
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Fig. 3. Density dependence of the neutrino emissivity Q 
in the NS core at T = 3 x 10 8 K and EOS A. Solid, dashed 
and dot-and-dashed lines lp, 2p, and 3p are plotted for 
proton superfluidity models displayed in Fig. 2. The dot- 
ted line shows Q in the presence of proton superfluidity 
lp and neutron superfluidity Int. The thick solid line is Q 
in non-supcrfluid matter. 
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Fig. 4. Density dependence of the critical temperatures 
of superfluidity of neutrons (models Ins, 2ns, and 3ns) in 
the NS crust (cf. Fig. 1). Vertical dotted line indicates 
neutron drip point. Horizontal line is the temperature of 
matter, T = 3 x 10 s K, adopted in Fig. 5. 
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Fig. 5. Density dependence of the neutrino emissivity Q 
in the NS crust at T = 3 x 10 8 K. Thin solid, dashed, 
and dot-and-dashed lines Ins, 2ns, and 3ns are calculated 
for neutron superfluidity models displayed in Fig. 4. The 
thick solid line is for non-superfluid matter. The core-crust 
boundary is artificially shifted to higher p for a better 
visualization of the effects of neutrino emission associated 
with 1 So Cooper pairing of neutrons. 

scopic theories which adopt a moderately strong medium 
polarization of proton-proton interaction. 

Finally, the dotted curve in Figs. |l| and | shows T cnt (p) 
for lnt m odel of 3 P2 neutron pairing (used by Kaminker et 
al. 2001 ). The curve falls within a wide scatter of T cn t(p), 



provided by microscopic theories. 

4. Main regulators of neutron-star cooling 

Theoretical cooling curves depend on many features of 
NS models. Below we discuss four main regulators of the 
cooling of middle-aged NSs: (1) central kernels, where neu- 
trino emission is enhanced by direct Urea process; (2) the 
effects of nucleon superfluidity on neutrino emission; (3) 
surface envelopes of light elements; and (4) surface mag- 
netic fields. The effects of these factors are illustrated in 
Figs. |-|7| The effects (1) and (2) are the strongest, al- 
though (3) and (4) can also be important. 

4-1- Direct Urea process in non-superfluid NSs 

The NS cooling is strongly affected by the presence of the 
central NS kernels (p > pr>, Table 1), where direct Urea 
process - the most powerful neutrino emission mechanism 
- is open. The effect is especially pronounced in the ab- 
sence of nucleon superfluidity. For instance, Fig. [| shows 



the density profile of the total neutrino emissivity through- 
out the NS core at a temperature T = 3 x 10 8 K. The 
emissivity of non-superfluid matter is shown by the thick 
line. The jump by about 7 orders of magnitude is asso- 
ciated with the direct Urea threshold, pr>. The NSs with 
the central density p c < po have no kernels with the en- 
hanced neutrino emission and show slow cooling. The NSs 
with p c > pd will have these kernels and show fast cool- 
ing. The presence of the kerne ls has dramatic effects on 
the NS cooling (Lattimer et al. 1991 ). 



4-2. The effects of superfluidity on neutrino emission 

First we discuss the effects of superfluidity in the NS cores. 
The adopted models of proton and neutron critical tem- 
peratures (Sect. 3) in the core are shown in Figs, [l] and 
0. The appropriate neutrino emissivities are displayed in 
Fig. 3. A proton superfluidity drastically changes the den- 
sity profiles of the neutrino emissivity, Q(p), as compared 
with the non-superfluid case (lines lp, 2p, and 3p in Figs. 
|and §. 

At not too high density, p ^ pi (the threshold den- 
sity pi is specified in Sect. 6), a strong proton super- 
fluidity almost switches off all neutrino processes involv- 
ing protons, particularly, modified Urea process and even 
the most powerful direct Urea process. In this regime, 
the main neutrino emission is provided by neutrino-pair 
bremsstrahlung in neutron-neutron scattering; its level is 
lower than in a slowly cooling non-superfluid NS regulated 
by non-suppressed modified Urea process. Since strong 
proton superfluidity 3p extends to higher p, this regime 
persists to higher p, for this superfluidity. 

At rather high p pn (the value pn is also specified in 
Sect. 6) the proton superfluidity dies out and we have the 
neutrino emissivity fully enhanced by direct Urea process. 
Again, pn is higher for proton superfluidity 3p. 

At pi < p < pn the neutrino emissivity is determined 
by direct Urea process partly suppressed by proton super- 
fluidity. Evidently, this regime is absent in non-superfluid 
NSs. In this regime, the emissivity is sensitive to the de- 
creasing slope of the T cp (p) profile and to EOS in the NS 
core. 

Thus, in any Q(p) profile, we can distinguish three 
distinctly different parts: (I) p ^ pi; (II) pi < p <, pn; 
(III) p <; pn. They will lead (Sect. 6) to three types of 
cooling superfluid NSs. 

Now we switch on neutron superfluidity in the core 
(model lnt) in addition to the proton one. For T = 3 x 
10 8 K, this superfluidity occurs in a density range from 



3 x 10 14 g cm~ 3 to 



6 x 10 14 g cm" 3 (Fig. 2). It 



creates a powerful splash of neutrino emission associated 
with Cooper pairing of neutrons as shown in Fig. 3 by the 
dotted line. This splash is actually the same for any model 
of proton superfluidity (lp, 2p, or 3p). 

In addition, we discuss the effects of the neutron su- 
perfluidity in the crust. The density profiles of the critical 
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Table 3. Surface temperatures of eight isolated middle-aged neutron stars inferred from observations 



Source 



RXJ0822-43 

IE 1207-52 

RXJ0002+62 

PSR 0833-45 (Vela) 

PSR 0656+14 

PSR 0633+1748 (Geminga) 

PSR 1055-52 

RX J1856-3754 



Igt 

[yr] 



3.57 

3.85 

3.95 b) 

4.4 C) 

5.00 

5.53 

5.73 

5.95 



[K] 



6.23Z 



6 10+ 005 

D - iu -0.06 



+0.03 



6.03 



5.83! 



5 96+ 02 

J - yD -0.03 



5.75 



+0.05 



j+0.03 
5 -0.04 



r 79 +0.05 
°- ' z -0.06 



Moder ) 



H 
H 
H 
H 

bb 
bb 
bb 



Confid. 
level 



95.5% 
90% 

95.5% 
68% 
90% 

90% 

d) 

d) 



References 



Zavlin et al. (1999) 



Zavlin et al. fll998| ) 
Zavlin & Pavlov ( |l99S| ) 



Pavlov et al. (2001) 



Possenti et al. (1996) 



Halpern & Wang (1997) 



Ogelman (1995) 



Pons et al. (2002) 



a ) Observations are interpreted either with a hydrogen atmosphere model (H), or with a black body spectrum (bb) 



^ The mean age taken according to Craig et al. (1997). 



c) According to Lyne et al. (|199q ) 
d > Confidence level is uncertain. 



Analytic fit with Si-ash atmosphere model of Pons et al. ( 2002 



temperature T cns (p) (models Ins, 2ns, and 3ns; Table 2; 
Figs. |l| and |^) and the appropriate neutrino emissivities 
Q{p) at T = 3 x 10 s K are shown in Fig. |. The thick 
line in Fig. exhibits the emissivity of the non-superfluid 
crust provided mainly by neutrino-pair bremsstrahlung in 
electron-nucleus scattering. One can observe strong ad- 
ditional peaks of neutrino emission from superfluid crust 
produced by Cooper pairing of neutrons. The peaks are 
very sensitive to the models of T cns (p). However, since the 
NS crust contains a small fraction of NS mass (about 1-5 
% or smaller, for the adopted NS models) this additional 
neutrino emission is important only in those NSs (Sect. 
6.3) whose cores produce low neutrino luminosity. 

4-3. Accreted surface layers 

A heat-blanketing envelope (Sect. 2) is traditionally as- 
sumed to be made of iron (e.g., Gudmundsson et al. 1983f ). 
However, the envelope may partly contain light elements 
(H, He) provided, for instance, by accretion. In the heat- 
blanketing layers, the heat is mainly carried by degener- 
ate electrons. The thermal conductivity of electrons which 
scatter off light ions is higher than the conductivity limited 
by scattering off iron ions. Thus, the accreted envelope is 
more heat transparent. The appropriate T s — Tb relation- 
ship was calculated by Potekhin et al. (1997). It depends 
on AM, the mass of accreted hydrogen or helium. The 
calculation takes into account nuclear burning of light ele- 
ments in a hot and/or dense plasma, i.e., AM is the total 
accreted mass which may partly be processed into heavier 
elements due to nuclear transformations in accreted mat- 



ter. The dependence of the photon NS luminosity on AM 
for four values of the temperature at the bottom of the 
heat-blanketing envelope (T h = 3 x 10 7 , 10 8 , 3 x 10 s , and 
10 9 K) is shown in Fig. 6. If the heat-blanketing envelope 
is fully replaced by accreted material, the NS luminosity 
can increase by about one order of magnitude, i.e., the 
surface temperature T s can increase by a factor of 2. 

4-4- Surface magnetic field 

The effects of strong magnetic field on the heat conduction 
through the blanketing NS envelope have been reviewed, 
for instance, by Ventura & Potekhin ( 2001 ). These effects 
depend on p, T, and the field strength B. One can distin- 
guish (i) classical effects of electron Larmor rotation, and 
(ii) the effects of quantization of electron orbits across the 
field lines (Landau states). 

The classical effects require lower magnetic fields. They 
become pronounced at w g r ^ 1, where u> g is the electron 
gyrofrequency, and r is the effective electron relaxation 
time. Then the conduction becomes anisotropic. The elec- 
tron thermal conductivity along the field lines remains 
equal to its non-magnetic value kq, while the conductivity 
in the perpendicular direction becomes kq/[1 + (w g r) 2 ]. 

The conditions of electron quantization are determined 
by characteristic density ps and temperature Tg, 



PB 

Tb 



7.045 x 10 3 (A/Z)B 3 1 i 2 gem" 

1.343 xl0 8 Si2 
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Fig. 6. Photon surface luminosity (redshifted as detected 
by a distant observer, left vertical axis) or redshifted effec- 
tive surface temperature (right vertical axis) of a "canon- 
ical" NS model (M = 1.4 Mq, R= 10 km) for four values 
of the "internal" temperature T, (at the bottom of the 
heat blanketing layer, pb = 10 10 g cm" 3 ) versus the mass 
of accreted material, AM. 



where Z is the mean charge number of atomic nuclei, A is 
the number of nucleons per a nucleus, pe = p/10 6 gem -3 , 
and B12 = B/10 12 G. The magnetic field remains non- 
quantizing as long as T > Tb- 

If p ii pb and T Tb, then the magnetic field acts 
as weakly quantizing. In this regime, the conductivities 
exhibit de Haas - van Alphen oscillations around their 
classical values. The field becomes strongly quantizing at 
p Pb and T <, Tb; it substantially modifies all compo- 
nents of the conductivity tensors well as the equation of 
state. 

The thermal structure of the NS envelopes composed 
of iron at B up to 10 16 G was studied, e.g., by Potekhin & 



Yakovlev fl200l[ ). Consider, for instance, the dipole field. 
Near the magnetic pole, the heat is carried away from the 
NS interior along the field lines. The electron-quantization 
effects amplify this conduction and make the polar regions 
of the heat-blanketing envelope more heat-transparent, in- 
creasing the local effective surface temperature for a given 
Tb. On the contrary, near the magnetic equator the heat 
propagates across the magnetic field. If w g r » 1, the 
classical Larmor-rotation effects strongly reduce the corre- 
sponding transverse conductivity. The equatorial regions 
become less heat transparent, which lowers the local ef- 
fective temperature for a given Tb. In an arbitrary ele- 
ment of the envelope, the heat conduction is affected by 



Fig. 7. Redshifted photon surface luminosity (left vertical 
axis) and mean effective temperature (right vertical axis) 
of a "canonical" NS with a dipole magnetic field, for two 
values of Tb and three models of the heat-blanketing en- 
velope (accreted mass AM = 0, lCr 12 M , and 1CT 7 M ) 
versus magnetic field strength B at the magnetic pole. 



both (longitudinal and transverse) conductivities. The NS 
photon luminosity is obtained by integration of the local 
radiated flux over the entire NS surface - see Potekhin & 



Yakovlev (2001) and references therein for the case of the 
envelopes composed of iron. Here, we supplement these 
results with new results obtained for partly accreted en- 
velopes with magnetic fields. The details of the new cal- 
culations will be presented elsewhere. Figure |?] displays 
the photon luminosity versus B for two selected values of 
Tb and three selected values of AM. The magnetic field 
affects the luminosity at B il 3 x 10 10 G. In the range 
of B up to B ~ 3 x 10 13 G, the equatorial decrease of 
the heat transport dominates, and the NS luminosity is 
lower than at B = 0. For higher B, the polar increase of 
the heat transport becomes more important, and the mag- 
netic field increases the photon luminosity. The variation 
of the photon luminosity by the magnetic field B ^ 10 15 
G does not exceed a factor of 30. 

The joint effect of the accreted envelope and the mag- 
netic field is demonstrated by the dot-dashed and dashed 
lines. As in the non-magnetic case, the accreted material 
makes the envelope more heat-transparent, thus increasing 
the luminosity at given Tb. Therefore, at B ~ 10 10 -10 13 
G, the magnetic field and the accreted envelope affect the 
thermal insulation in opposite directions. At higher field 
strengths, both effects increase the luminosity. However, 
as evident from Fig. 0, the dependence of this increase on 
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B and AM is complicated. In particular, at B <; 10 14 G, 
the effect of the accreted envelope is generally weaker than 
in the non-magnetic case. 

5. Observational data 

We will confront theoretical cooling curves with the results 
of observations of thermal emission from eight middle- 
aged isolated NSs. In discussing the observations and their 
theoretical interpretation we closely follow the considera- 
tion in KYG. 

The observational data are the same as KYG. They are 
summarized in Table 3 and displayed in Figs. 8-13. The 
three youngest objects (RX J0822-43, IE 1207-52, and 
RX J0002+62) are radio-quiet NSs in supernova remnants. 
The oldest object, RX J1856-3754, is also a radio-quiet 
NS. The other objects, Vela, PSR 0656+14, Geminga, and 
PSR 1055-52, are observed as radio pulsars. The NS ages 
are either pulsar spindown ages or the estimated super- 
nova ages. The age of RX J1856-3754 was estimated by 
Walter (2001) from the kinematical data (by identifying 
a possible companion in the binary system existed before 
the supernova explosion). We use the value t — 9 x 10 5 
yr mentioned in the subsequent publication by Pons et al. 
( p002l ). 

For the four youngest sources, the effective surface 
temperatures T s °° are obtained from the observed X-ray 
spectra using hydrogen atmosphere models. Such mod- 
els are more consistent (e.g., Pavlov & Zavlin 2002 ) with 
other information on these sources (distances, hydrogen 
column densities, inferred NS radii, etc.) than the black- 
body model of NS emission. On the contrary, for the next 
three sources we present the values of T s °° inferred using 
the blackbody spectrum because the blackbody model is 
more consistent for these sources. Finally, for RX J1856- 
3754 we adopt the values inferred using t he ana lytic fit 
with Si-ash atmosphere model of Pons et al. ( [2002 ). We ex- 
pect that the large error bar of T s °° provided by this model 
reflects poor understanding of thermal emission from this 
source (e.g., Pons et al. 2002, Burwitz et al. 2001, Gansicke 
et al. [200lj Kaplan et al. |200l| , Pavlov et al. |2002D . 



6. Theory and observations 

6.1. N on- superfluid neutron stars 

There is a large scatter of observational limits on T s °° for 
the eight sources. Three sources, the youngest RX J0822- 
43, and two oldest, PSR 1055-52 and RX J1856-3754, 
seem to be hot for their ages, while the other ones, espe- 
cially Vela and Geminga, look much colder. We can inter- 
pret all observational data with the cooling curves using 
the fixed (the same) EOS and models of the critical tem- 
peratures T c (p) (Sect. 3) for all objects. The results are 
presented in Figs. || |l3|. 

In the absence of any superfluidity there are two well- 
known, distinctly different cooling regimes, slow and fast 



cooling. The slowly cooling NSs have low masses, M < 
Md . They mainly loose their energy via neutrino emission 
in modified Urea processes. For given EOSs in the NS core, 
the cooling curves of middle-aged NSs are almost the same 
for all M from about M Q to Md (e.g., Page & Applegate 



1992, Gnedin et al. 2001), being not very sensitive to EOS. 



The fast cooling occurs via a very powerful direct Urea 
process (Lattimer et al. |1991| ) if M > M D + 0.003 M Q . 
The cooling curves are again not too sensitive to the mass 
and EOS. The middle-aged rapidly cooling NSs are much 
colder than the slowly cooling ones. Two examples, 1.35 
and 1.5 M Q non-superfluid NSs (EOS A), are displayed in 
Fig. ||. The transition from the slow to fast cooling takes 
place in a very narrow mass range, Md < M Md+ 
0.003 Mq. In order to explain these observational limits 
with the non-superfluid NS models one should make an 
unlikely assumption that at least the masses of the Vela 
and Geminga pulsars fall in that narrow mass range. Thus, 
several sources exhibit the intermediate case between the 
slow and fast cooling. In the absence of superfluidity, this 
is highly unlikely. 

6.2. Proton superfluidity and the three types of cooling 
neutron stars 

The next step is to interpret the observations (Fig. p|-p^|) 
by cooling of superfluid NSs. It turns out that various su- 
pcrfluids affect NS cooling in different ways. Our main as- 
sumptions would be that the proton superfluidity is rather 
strong at p <, p^>, while the 3 P2 neutron superfluidity is 
rather weak (Sect. 6.4). We start with the effects of proton 
superfluidity. 

The dot-and-dashed cooling curves in Figs. ||-|ll] are 
computed assuming the proton superfluidity alone. The 
proton pairing lp is adopted in Fig. 2p in Fig. [To|, and 
3p in Fig. O. EOS A is used in the models in Figs. |9j and 
[U| and EOS B in Fig. [HI 

As seen from Figs. ?rpl the proton superfluidity 
leads to the three cooling regimes (instead of two in non- 
superfluid NSs): slow, moderate, and fast. They reflect 
three neutrino emission regimes from the NS cores dis- 
cussed in Sect. 4.2 (Fig. |3|). Accordingly, there are three 
types of cooling NSs with distinctly different properties. 

(I) Low-mass, slowly cooling NSs. Their central densi- 
ties p c and masses M = M(p c ) obey the inequalities 



A 



< 



Pi- 



M <, Mi = M(pi) 



(9) 



where p\, introduced in Sect. 4.2, is discussed below in 
more detail. 

(II) Medium-mass, moderately cooling NSs, with 



Pi SS Pc £ pa, Mi <S M < M„ = M(p„), 



(10) 



where pu, also introduced in Sect. 4.2, is specified below. 
(Ill) Massive, rapidly cooling NSs, 



pu & Pc < 



Mu ^ M < M„ 



(11) 
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Fig. 8. Observational limits on surface temperatures of 
eight NSs (Table 3) compared with cooling curves of non- 
superfluid NSs (EOS A) of M = 1.35 and 1.5 M . 



where p cm ax and M max refer to the maximum-mass con- 
figuration (Table 1). 

The threshold values of piji and Mj t n depend on a 
proton superfluidity model T cp (p), EOS in the NS core, 
and on a NS age. The main properties of regimes I, II, 
and III are as follows. 

(I) The slowly cooling NSs are those where the direct 
Urea process is either forbidden by momentum conserva- 
tion (p c < pd) or greatly suppressed by the strong proton 
superfluidity (Fig. ||). 

In particular, the cooling is slow for p c < pn and 
M < Mrj in the absence of proton superfluidity. This is 
the well-known ordinary slow cooling which is mainly reg- 
ulated by the neutrino emission produced by the modified 
Urea process (Sect. 6.1). However, for the conditions dis- 
played in Figs. the proton superfluidity in regime 
I is so strong that it almost switches off both, the modi- 
fied Urea process everywhere in the NS core and the di- 
rect Urea process at p > pr>- Then the neutrino emis- 
sion is mainly produced by neutrino bremsstrahlung in 
neutron-neutron scattering (Sect. 4.2) unaffected by the 
neutron superfluidity in the NS core (that is assumed to 
be weak). The bremsstrahlung is less efficient than the 
modified Urea process and leads to a slower cooling than 
in a non-supcrfluid NS. It can be referred to as the very 
slow cooling. 

The analysis shows (KYG) that, for the selected mod- 
els, the regime of very slow cooling holds as long as the 




Fig. 9. Observational limits on surface temperatures of 
eight NSs compared with cooling curves of NSs (EOS A) 
with masses from 1.35 to 1.55 Mq (after KYG). Dot-and- 
dashed curves are obtained including proton superfluidity 
lp alone. Solid curves include, in addition, model Ins of 
neutron superfluidity. Dotted lines also take into account 
the effect of neutron superfluidity Int. 

proton critical temperature in the NS center is higher than 
a threshold value: 



t cp ( Pc )Zt£\ Pc ). 



(12) 



Comparing the neutrino emissivities of the indicated reac- 
tions (e.g., Yakovlev et al. 1999), KYG obtained the simple 



estimates: T c ( p ; (p) - 5.5 T for p < p D , and T^'(p) - 17 T 
for any p which is several per cent higher than pn, where 
T is the internal NS temperature. There is a continuous 
transition (Fig. 9 in KYG) between the presented values 
of T c ( p' (p) in the narrow density range near prj. 

To make the analysis less abstract we notice that in 
a very slowly cooling NS the internal temperature is T ~ 
5.5 x 10 s , 4 x 10 8 , and 1.5 x 10 s K at t ~ 4 x 10 3 , 2.5 x 10 4 , 
and 4 x 10 5 yr, respectively. 

Now we can explicitly specify the maximum central 
densities p\ and masses Mj of very slowly cooling NSs in 
Eq. (g) for the cases of study, in which p\ > pn (Mi J> 
Md). If p c > Pd, then we can define p\ as the density 
value on the decreasing, high-density slope of T cp (p) (Figs. 

|l| and||) which corresponds to T cp (pj) = Top (pi). It gives 
the central density of the NS and M\ = M(p{). 

For the conditions displayed in Figs. p| ]IT|, the cooling 
curves (dot-and-dashed lines) of all low-mass NSs with any 
proton superfluidity (lp, 2p, or 3p) are very similar. For 
instance, the 1.35 M Q curve in Fig. || is plotted just as an 



(i), 
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Fig. 10. Observational limits on surface temperatures of 
NSs compared with cooling curves of NSs (EOS A) with 
several masses M in the presence of proton superfluidity 
2p (from KYG). Dot-and-dashed curves are obtained as- 
suming non-superfluid neutrons. Solid curves include, in 
addition, model Ins of neutron superfluidity. 3 P2 neutron 
pairing is neglected. 
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Fig. 11. Observational limits on surface temperatures of 
NSs compared with cooling curves of NSs (EOS B) with 
masses from 1.3 to 1.6 M© (from KYG). Dot-and-dashed 
curves are obtained using model 3p of proton superfluidity. 
Solid curves include, in addition, model Ins of neutron 
superfluidity. 



example; all the curves are almost identical in the mass 
range M Q ^ M < M\. Moreover, the cooling curves are 
not too sensitive to EOS and are insensitive to the exact 
values of the proton critical temperature T cp , as long as the 
inequality ([l2]) holds. These curves are noticeably higher 
than the curves which describe the ordinary slow cooling 
in the absence of superfluidity (e.g., Fig. |J). It is remark- 
able that they seem to be almost model-independent. 

For the conditions in Figs. p|-pd|, the three relatively 
hot sources, RX J0822-43, PSR 1055-52, and RX J1856- 
3754, can be treated as these very-slow-cooling (low-mass) 
models. 

(II) The moderately cooling stars are the NSs which 
possess central kernels where the direct Urea process is 
allowed but moderately suppressed by proton superfluid- 
ity. The existence of a representative class of these NSs is 
solely due to proton superfluidity. 

According to KYG, for the selected cooling models, 
the proton critical temperature in the center of a medium- 
mass NS should roughly satisfy the inequality 



Tg l \p c )£T cp (p c )<Tg\p c ), 



(13) 



with T, 



(ID 



cp 



3 T . Thus, one can find pu which corresponds 
to T cp (pu) — Tc^^u) and determines Mn, the maximum 
mass of moderately cooling NSs in Eq. ( [Toj ) . 



The surface temperatures of these medium-mass NSs 
are governed by proton superfluidity in the NS central 
kernels, p > p\. One can observe (Figs. (9pf[l]) a steady 
decrease of the surface temperatures with increasing M. 
Fixing the proton superfluidity and EOS, one can deter- 



mine (Kaminker et al. 2001 and Yakovlev et al. 2001b) the 



masses of moderately cooling NS, which means "weigh- 
ing" NSs. In this fashion one can weigh five isolated NSs 
(IE 1207-52, RX J0002+62, Vela, PSR 0656+14, and 
Geminga) as shown in Figs. For instance, adopt- 

ing EOS A and proton superfluid lp (Fig. ||) KYG obtain 
the masses in the range from ss 1.36 M (for IE 1207-52) 
to « 1.465 M Q (for Vela and Geminga). For EOS A and 
proton superfluid 2p (Fig. [Io|) KYG obtain higher masses 
of the same sources. Obviously, the properties of moder- 
ately cooling NSs are extremely sensitive to the decreasing 
slope of T cp (p) in the density range from pi to pu (and in- 
sensitive to the details of T cp (p) outside this range). 

(Ill) Massive NSs show fast cooling similar to the fast 
cooling of non-superfluid NSs. These stars have central 
kernels where the direct Urea process is either unaffected 
or weakly suppressed by the proton superfluidity. In such 



kernels, T cp (p) 



< T t 



(ii) 



cp 



The central densities and masses 



of these NSs lie in the range given by Eq. (|1 1|) . Their 
thermal evolution is not very sensitive to the model of 
T cp (p) and to EOS in the stellar core. Note that if /? C max < 
pu, the rapidly cooling NSs do not exist. In the frame of 
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this interpretation, no NS observed so far can be assigned 
to this class. 



6.3. Crustal superfluidity and slow cooling 

At the next step, we retain proton superfluidity and add 
neutron superfluidity Ins in the NS crust. These mod- 
els are shown by the solid curves in Figs. §-0. For the 
moderately or rapidly cooling middle-aged NSs, they are 
fairly close to the dot-and-dashed curves. This is quite ex- 
pected (e.g., Gnedin et al. |200l| ): the 1 So neutron superflu- 
idity is mainly located in the NS crust which is much less 
massive than the NS core. Thus, the crustal superfluid- 
ity does not affect noticeably the proposed interpretation 
of IE 1207-43, RX J0002+62, Vela, PSR 0656+14, and 
Geminga in terms of moderately cooling NSs. 

However, as pointed out by Yakovlev et al. ( 2001b| ), 
this crustal superfluidity strongly affects the slow cool- 
ing of low-mass NSs, and the effects are twofold. First, at 
t 5s 3 x 10 5 yr the neutrino luminosity due to 1 Sq pair- 



ing of neutrons may dominate the sufficiently low neutrino 
luminosity of the stellar core. Second, at t <; 10 5 yr the 
1 So neutron superfluidity reduces the heat capacity of the 
crust. Both effects accelerate NS cooling and decrease T s °° 
(Figs. p|-pd|) violating the interpretation of the three suffi- 
ciently hot sources, RX J0822-43, PSR 1055-52, and RX 
J1856-3754. The interpretation of RX J1856-3754 is af- 
fected to a lesser extent, as a consequence of the rather 
large error bar of T s °° for this source (Sect. 5). 

However, the interpretation can be rescued by the ap- 
propriate choice of T cns (p). It is sufficient to focus on the 
interpretation of RX J0822-43, PSR 1055-52, and RX 
J1856-3754, as the very slowly cooling NSs. For certainty, 
we take EOS B, M = 1.3 M Q , and proton superfluid 3p. 
The results are presented in Figs, [l^ and p]| The dot- 
and-dashed line is the same as in Fig. [ll] and neglects 
the crustal neutron superfluidity. Thick solid line is also 
the same as in Fig. O. It includes an additional effect of 
crustal superfluid Ins and lies below the observational lim- 
its on T s °° for the sources in question (or almost below in 
case of RX J1856-3754). To keep the proposed interpreta- 
tion of the three sources one must raise the cooling curves 
calculated including the crustal superfluidity. To this aim, 
one should suppress the neutrino emission associated with 
1 So pairing of neutrons (Fig. ||). Recall that in a middle- 
aged NS this emission is mainly generated (Fig. ||) in two 
relatively narrow layers, near the neutron drip point and 
near the crust-boundary interface, where the local NS tem- 
perature T is just below T cns (p). Since the Cooper-pairing 
neutrino luminosity is roughly proportional to the widths 
of these emitting layers, one can reduce the luminosity by 
reducing the widths. This can be done by setting T™g X 
higher and by making the sharper decrease of T cns (p) in 
the wings. 

For example, taking crustal superfluid 2ns instead of 
Ins (Figs. | and |) one obtains the dashed cooling curve 



in Fig. |l2| which comes much closer to the dot-and-dashed 
curve than the thick solid curve (model Ins). Note that 
the cooling curves are insensitive to the details of T cns (p) 
profile near the maximum, as long as T™ x ^ 5 x 10 9 K, 
but they are extremely sensitive to the decreasing slopes of 
Tcns(p)- On the other hand, by taking the smoother and 
lower T cns (p), model 3ns, one obtains a colder NS than 
needed for the interpretation of the observations (long- 
dash line in Fig. |l^). Therefore, 1 So neutron superfluidity 
with maximum T™g X < 5 x 10 9 K and/or with smoothly 
decreasing slopes of the T cns (p) profile near the crust-core 
interface and the neutron drip point violates the proposed 
interpretation of the observational data. 

Moreover, the observations of RX J0822-43, PSR 
1055-52, and RX J1856-3754 can be fitted even with the 
initial model Ins of the crustal superfluidity (Fig. |l3| ). 
The high surface temperature of RX J0822-43 can be ex- 
plained assuming additionally the presence of a low-mass 
(2 x 10 -11 Mq) heat-blanketing surface envelope of hydro- 
gen or helium. This effect is modeled using the results of 
Potekhin et al. (|1997[) described in Sect. 4 (Fig. §). Light 



elements raise T s °° at the neutrino cooling stage (curve 
ace in Fig. |l3| ). In order to explain the observations of 
PSR 1055-52 and RX J1856-3754, one can assume again 
model Ins of crustal superfluidity, iron surface and the 
dipole surface magnetic field (~ 10 12 G at the magnetic 
pole; line mag in Fig. [l3]). Such a field makes the NS sur- 
face layers overall less heat-transparent (Fig. Q), rising T s °° 
at t ^ 3 x 10 5 yr. Note that the dipole field k, 3 x 10 13 G 
has the opposite effect, resembling the effect of the surface 
envelope of light elements. 

To summarize, one can additionally vary cooling 
curves by assuming the presence of light elements and the 
magnetic field on the NS surface (line acc-mag in Fig. |l3| ). 
However, these variations are less pronounced than those 
due to nucleon superfluidity. For instance, one cannot rec- 
oncile the cooling curves with the present observations of 
PSR 1055-52 assuming model 3ns of the crustal superflu- 
idity with any surface magnetic field. 

6.4- 3 P2 pairing of neutrons in the NS core 

Now we focus on the 3 P2 neutron pairing neglected so far. 
Its effects are illustrated in Fig. [| as an example. They are 
qualitatively similar for the other cooling models in Figs. 
and In Fig. § we take the cooling models obtained 
including proton superfluidity lp and crustal superfluid- 
ity Ins, and add the 3 P2 neutron superfluidity (model lnt, 
Table 2) in the core. This gives the same (solid) cooling 
curves for the young NSs which have the internal temper- 
atures T above the maximum value of T™^ x « 3 x 10 s K. 
However, when T falls below T™j* x , one obtains (dots) a 
strong acceleration of the cooling associated with the pow- 
erful neutrino emission due to 3 P2 neutron pairing (Fig. 
H). This emission greatly complicates the proposed inter- 
pretation of older sources, PSR 0656+14, Geminga, PSR 
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Fig. 12. Cooling curves of 1.3 M Q NS model (EOS B) 
versus observations of RX J0822-43, PSR 1055-52, and 
RX J1856-3754 (after KYG). Dot-and-dashed curve: pro- 
ton superfluidity 3p in the NS core. Solid, short-dashed, 
and long-dashed curves include, in addition, models Ins, 
2ns, and 3ns of crustal neutron superfluidity, respectively. 
Thick solid line is the same as in Fig. 01 



1055-52, and RX J1856-3754. The complication arises for 
a wide class of T CDt (p) profiles with T™f x from ~ 10 8 K to 
~ 3 x 10 9 K (in this way, it can be regarded as model- 
independent). The Cooper-pairing neutrino emission of 
neutrons induces really fast cooling of older sources even 
if their mass is low, M < Md (Sect. 6.2). To avoid this 
difficulty one can assume (Kaminker et al. 2001) weak 3 P2 
pairing, T cnt (p), with maximum T™f x < 10 8 K; it does 
not affect the proposed interpretation. 



7. Summary 

Following KYG we summarize the effects superfluids on 
NS cooling. 

(a) Strong proton superfluidity in the NS cores, com- 
bined with the direct Urea process at p > pu, separates 
(Sect. 6.2) the cooling models into three distinctly differ- 
ent types: (I) slowly cooling, low-mass NSs (M ^ Mi); (II) 
moderately cooling, medium-mass NSs (Mi <, M <, Mn); 
(III) rapidly cooling, massive NSs (M J; -Mn). The regime 
of moderate cooling cannot be realized without the proton 
superfluidity. 

(b) Strong proton superfluidity in the NS core is re- 
quired to interpret the observational data on the three 
sources, RX J0822-43, PSR 1055-52, and RX J1856-3754, 
hot for their ages, as the very slowly cooling NSs (Sects. 



Fig. 13. Cooling curves of 1.3 M Q NS model versus ob- 
servations of RX J0822-43, PSR 1055-52, and RX J1856- 
3754. Dot-and-dashed curve: proton superfluidity 3p in the 
NS core. Solid curves include, in addition, model Ins of 
crustal neutron superfluidity. Thick solid line is the same 
as in Figs. |ll| and |lj. Thin solid curve ace is calculated 
assuming the presence of 2 x 10 _11 M Q of hydrogen on 
the NS surface. Thin solid curve mag is obtained assum- 
ing the iron surface and dipole magnetic field (10 12 G at 
the magnetic pole). Thin solid curve acc-mag is obtained 
assuming both, the magnetic field (10 12 G) and accreted 
envelope (2 x 10" n M Q ). 



6.2 and 6.3). Within this interpretation, all three sources 
may have masses from about M Q to Mi; one cannot de- 
termine their masses exactly or distinguish EOS in the NS 
core from the cooling models. 

(c) Strong proton superfluidity is needed to inter- 
pret observations of the other sources, IE 1207-52, RX 
J0002+62, Vela, PSR 0656+14, and Geminga, as the 
medium- mass NSs. One can "weigh" these NSs, i.e., de- 
termine their masses, for a given model of T cp (p) and a 
given EOS. The weighing is very sensitive to the decreas- 
ing slope of T cp (p) in the density range pi Si p ^ pu (Sect. 
6.2). 

(d) Strong or moderate 3 P2 neutron superfluidity in 
the NS core initiates rapid cooling due to the neutrino 
emission resulted from neutron pairing. This would inval- 
idate the proposed interpretation of the old sources like 
PSR 0656+14, Geminga, PSR 1055-52, and RX J1856- 
3754. One can save the interpretation assuming a weak 
3 P 2 neutron superfluidity, T c ™ ax < 10 8 K (Sect. 6.4). 

(e) 1 So neutron superfluidity in the crust may initi- 
ate a strong Cooper-pairing neutrino emission, decrease 
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substantially T s °° of the slowly cooling NSs, and weaken 
the interpretation of RX J0822-43, PSR 1055-52, and RX 
J1856-3754 (although it does not affect significantly the 
moderate or fast cooling). The interpretation can be saved 
by assuming that the maximum of T cns (p) is not too small 
(T™ x £ 5 x 10 9 K) and the profile of T cns (p) decreases 
sharply in the wings (Sect. 6.3). 

(f) The interpretation of the slowly cooling sources is 
sensitive to the presence of the surface magnetic fields 
and/or heat-blanketing surface layer composed of light el- 
ements (Sect. 6.3). 

(g) No isolated middle-aged NSs observed so far can 
be idenfied as a rapidly cooling NS. In the frame of the 
proposed models, these NSs do not exist for those EOSs 
and T cp (p) profiles for which M max < M\\. 

If the proposed interpretation is correct, one can make 
the following conclusions on the properties of dense matter 
in NS interiors. 

(i) Strong proton superfluidity is in favor of a not 
too large symmetry energy at supranuclear densities 
(Kaminker et al. 2001). A very large symmetry energy 



would mean a high proton fraction which would suppress 
proton pairing. On the other hand, the symmetry energy 
should be not too small to open the direct Urea process 
at p > p D . 

(ii) Weak 3 P2 neutron pairing is in favor of a not too 
soft EOS in the NS core (Kaminker et al. |200l| ). The soft- 
ness would mean a strong attractive neutron-neutron in- 
teraction and, therefore, strong neutron pairing. 

(iii) The adopted features of the crustal neutron super- 
fluidity are in favor of those microscopic theories which 
predict T cna (p) profiles with T^f > 5 x 10 9 K. This is 
in line with many microscopic calculations of the super- 
fluid gaps which include the medium polarization effects 
in neutron- neutron interaction (e.g., Lombardo & Schulze 
2001). However, the reduction of the gap by the medium 



polarization should not be too strong, and the decreasing 
slope of T cns (p) should be rather sharp. These require- 
ments constrain the microscopic theories. 

The proposed interpretation of the observations relates 
the inferred NS masses to the superfluid properties of NS 
interiors. By varying EOS and the proton critical tem- 
perature, one can attribute different masses to the same 
sources. If, on the other hand, one knew the range of 
masses of the cooling middle-aged NSs it would be possi- 
ble to draw definite conclusions on the superfluid state of 
their interiors, first of all, on the proton critical tempera- 
ture, T'cp(p)- 

The presented analysis may seem too simplified be- 
cause it neglects a possible presence of other particles in 
the NS cores (muons, hyperons, quarks). It is expected 
that the inclusion of other particles and the effects of su- 
perfluidity of hyperons or quarks will complicate theoret- 
ical analysis but will not change the basic conclusion on 
the existence of the slowly, moderately, and rapidly cool- 
ing NSs. 



The calculations show that the cooling of middle-aged 
NSs with M < M\ is sensitive to the density profile of free 
neutrons near the crust bottom and neutron drip point. 
The presented calculations used only one model of the 
free- neutron distribution in the crust, assuming spherical 
atomic nuclei at the crust bottom. It would be interesting 
to consider the models of crust matter with non-spherical 
nuclei (e.g., Pethick & Ravenhall 1995) and the effects of 



superfluidity of nucleons confined in the atomic nuclei in 
the NS crust. 

The determination of T s °° from observational data is a 
very c omplic ated problem (as described, e.g., by Pavlov & 
Zavlin pOOl ). It requires very high-quality data and theo- 
retical models of NS atmospheres. Thus, the current values 
of T s °° may change substantially after the forthcoming ob- 
servations and new theoretical modeling. These changes 
may affect the proposed interpretation of the observa- 
tional data, first of all, of RX J0822-43, PSR 1055-52, 
and RX J1856-3754. For instance, RX J1856-3754 may 
have a colder surface (T s °° ~ 0.25 MK), than assumed in 



the above analysis, with a hot spot (e.g., Pons et al. 2002, 
Burwitz et al. 2001, Gansicke et al. [2001 ). If confirmed, the 



lower T s °° might be explained by the effect of P2 neutron 
pairing (Fig. |). 

The future observations of the thermal emission from 
these sources will be crucial for understanding the super- 
fluid properties of NS matter. 

The proposed theory can explain the existence of NSs 
within a broad range of T s °°. However, it would be unable 
to explain too hot and too cold objects (with T s °° essen- 
tially higher than the highest cooling curve and essentially 
lower than the lowest curve in Figs. Ull]). A discovery of 
such hot or cold isolated NSs would be of special interest. 
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